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ABSTRACT 


The  elastic  scattering  of  5.0  MeV  deuterons  incident  on  40Ca 
was  measured  for  the  forward  angles.  The  purpose  of  the  measurements 
was  to  test  the  sensitivity  of  the  optical  model  parameters,  particularly 
to  small  deviations  from  Rutherford  scattering.  The  problems  of  separat- 
ing  the  other  predominant  elastic  peaks,  0  and  C  contaminants  from 
40Ca  was  encountered.  Large  deviations  from  the  Rutherford  scattering 
cross-section  for  both  these  nuclei  were  attributed  to  an  appreciable 

nuclear  part  of  the  scattering  amplitude  predicted  by  the  optical  model. 
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CHAPTER  1 


INTRODUCTION 

The  purpose  of  the  thesis  is  to  report  measurements  of  the 
absolute  cross-sections  in  the  forward-angle  region  between  7  to  30 
degrees  for  deuterons  scattered  from  calcium  at  5.0  MeV  and  to  discuss 
the  results  of  the  oscillatory  deviations  predicted  in  the  optical-model 
calculations.  The  investigation  was  carried  out  because  of  the  interest 
in  the  deuteron-nucleus  interaction,  particularly  since  little  experi¬ 
mental  information  is  available  for  the  forward  angles.  The  thesis 
attempts  to  resolve  some  of  the  difficulties  and  suggest  comments  or 
solutions  for  parts  of  the  experiment  not  completed. 

Angular  distributions  of  5.0  MeV  deuterons  scattered  from  C, 

160  and  40Ca  nuclei  were  measured  between  7  to  160  degrees.  At  low  beam 

energies,  say  5.0  MeV,  one  would  expect  scattering  from  the  above  nuclei  to 

be  pure  Coulomb  in  the  forward  angles .  However ,  the  forward-scattering  ampli¬ 
tude  contains  not  only  a  Coulomb  part  but  also  nuclear  contributions. 

As  the  scattering  angle  approaches  zero  degrees,  the  centre  of  mass  of 
the  elastic  peaks  shift  with  angle:  12C  shifts  into  160  while  both  peaks 
move  into  40Ca  [Fig.  l.l].  One  cannot  subtract  the  12C  and  160  peaks 
from  40Ca  without  taking  into  account  the  appreciable  nuclear  part. 

For  incident- pro ton  energies  of  2.5,  3.0  and  3.5  MeV,  angular 
distributions  were  obtained  for  ratios  between  relative  and 
absolute  cross  sections.  Unfortunately,  the  cross  sections  for  the 
elastic  scattering  of  low-energy  protons  from  C  and  0  nuclei  deviate 
from  Rutherford  scattering. 
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Figure  1.1  Elastic  Peaks  at  17,  20,  24  and  30  degrees  for 
5.0  MeV  deuterons. 


CHAPTER  2 


THEORY 

2.1  Rutherford  Scattering 

When  incident  particles  are  scattered  from  a  target  nucleus, 
the  interaction  may  be  described  in  terms  of  their  charges.  Unless  the 
incident  particles  have  sufficient  energy  to  overcome  the  repulsive 
Coulomb  force  of  the  nucleus,  they  cannot  approach  closely  enough  to 

the  target  nucleus  to  be  affected  by  the  short-range  nuclear  force.  The 
nonrelativis tic  Rutherford  scattering  cross  section  is 

2 

dCT(0)  1  ,  p  t:  N  4  /Q/,N  ,  ,, 

d'n(e',>7  "  4  (-7F~)  csc  (0/2)  mb/sr  (2.1-1) 

The  scattering  angle  is  0  while  E  is  the  incident  laboratory  energy. 

The  charges  z^e  and  z^_e  are  of  the  particle  and  target  respectively. 
Since  the  scattering  into  <{)  is  isotropic 

Jq2TT  d(J>  =  2 IT  (2.1-2) 

then  the  solid  angle  reduces  to 


dft(0,<J))  =  2tt  sin0  d0 


(2.1-3) 


The  above  expression  increases  rapidly  for  small  angle  0  and  becomes 
infinite  as  0  approaches  zero  degrees.  Rewriting  in  terms  of  the  impact 


parameter  b. 


z  z  e 

db  =  -  (— — )csc2(e/2)  d9 


(2.1-4) 


shows  that  small-angle  scattering  corresponds  to  very  large  impact 
parameters.  Therefore,  the  nuclear  contributions  should  get  smaller  as 
the  impact  parameter  b  becomes  larger  and  the  Coulomb  force  predominates. 
Angular  distributions  may  be  plotted  as  the  ratio  (da(0) /d£3)  ^  /  (da(0) /d^)^ut^ 

which  effectively  divide  out  the  Coulomb  effects  and  show  graphically 
only  the  nuclear  part. 


3 


.  •  '  • 


4 


2.2  Optical  Model 

From  charged  particle  data,  the  optical  model  (H071)  is  a  use¬ 
ful  description  of  the  elastic-scattering  results,  and  the  aim  is  to  find 
smooth  trends  in  the  potential  parameters  as  a  function  of  mass  number 
and  energy.  The  optical  model  used  is  defined  by  the  sum  of  the  follow¬ 
ing  potentials: 

i)  Coulomb  potential 

gPV2  [3-(r/Rc>2!  for  r  <  R 

2R  C 

c 

ZpZ  e2  r 

t  for  r  >  R 

-  .  -  c 

r 

R  =  r  A1/3 
c  c 


The  Coulomb  potential,  would  be  that  produced  by  a  uniform 

charge  distribution  of  classical  radius,  Rc.  Since  the  experimental 

results  are  not  sensitive  to  the  value  of  R^,  throughout  the  analysis, 

R  is  fixed  at  R  =  1.3  A1/3  where  A  is  the  mass  number  of  the  target 
c  c 

nucleus . 

ii)  Central  potential, 

real  part  -Vf (x  ) 

imaginary  part  -W  f (x  )  -  4a  W_df(x  )/dr 

s  w  wD  w 

The  function  f(x^)  is  the  usual  Woods-Saxon  form  factor, 

f(x.)  =  [l  +  exp(x.)]  1  where  x.  =  (r-r . A1 / 3) /a.  with  the  appropriate 
i  l  i  l  i 

radius  and  dif fuseness geometrical  parameters,  (r  ,  a  ,  r  ,  a  ,  r 

o  o  w  w  so, 

a  ).  The  imaginary  potential  is  composed  of  a  volume  part  W  ,  (set 
so  s 


and 


* 
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equal  to  zero  in  the  present  work)  and  a  surface  part  W^#  The 
factor  4  appearing  is  required  so  that  the  surface  form  factor 
is  a  maximum  when 

4a  df (x  ) / dr/ x  =  1 
w  w 

iii)  Spin-orbit  potential 

(fi/m  c)V  (1/r)  df (x  )/dr  L*a 

7T  SO  SO  '*>  ~ 

where  Ch/m^c)2  =  2  fm2.  The  operator  rr  is  defined  in 

terms  of  the  spin  angular  momentum  S  as  follows: 

S  =  ^  (h/ 2 )q  proton 

ha  deuteron 

Finally,  r-r1  is  the  relative  separation  of  the  interacting  particles 
describing  a  non-local  potential  U(r,r')  represented  by  the  expression 
U(r-r’)  =  U[(r-r’ )/2]exp[-C(r-r' )/3}2]/tt3/2  g3 

Since  3  is  a  measure  of  the  range  of  non-locality,  in  the  low  energy 
region,  3^0  then  U(r,r')  U(r)  becomes  a  local  potential  (H071). 

*v  *v  ^  * 

The  above  non-local  potential  was  calculated  from  a  nucleon-nuclear 
interaction.  The  optical  model  takes  into  account  all  interactions 
between  the  incident  particle  and  the  target  nucleons. 

In  this  work  we  consider  the  elastic  scattering  of  incident 
particles  from  a  target  nucleus.  In  the  local  approximation,  the  momentum 

transferred  between  the  collisions  with  the  target  nucleons  must  be 

small  in  order  that  the  internal  state  of  the  nucleus  remains  unchanged. 

This  implies  that  long  range  collisions  or  Coulomb  scattering  may  contri¬ 
bute  strongly  to  forward  angle  scattering. 


. 


There  are  certain  ambiguities  in  the  optical  parameters 
occurring  in  the  analysis  of  the  experimental  data.  The  procedure 
is  to  begin  with  acceptable  parameters  and  then  iterate  to  optimize 
a  fit  to  the  measured  data.  The  iteration  is  carried  out  until 
the  discrepancy  between  the  experimental  and  the  theoretical  cross- 
sections  is  a  minimum.  However,  there  is  the  danger  of  obtaining 
an  unphysical  solution  for  the  above  method.  Such  ambiguities  are 
most  generally  associated  with: 

a)  the  incorrect  potential  form  factor; 

b)  Vr11,  where  n=l,  2.  or  3  or  Wa  =  const.; 

w 

c)  many  potential  depths  give  comparable 
fits  to  the  data;  and 

d)  over-all  normalization  error  for  the 


experimental  data. 


' 
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2.3  Scattering  Amplitudes 


For  a  nucleon-nucleus  interaction,  the  incoming  and  outgoing 
particles  of  spin  components  m  and  m1  satisfy  the  Schrodinger  equation. 


=  V 

and  the  solution  (Ho  71)  may  be /written  asymptotically  as 

H'm  r  5.  „  [1  +  12/ik(r-z)] 


(2.3-1) 


m 

s 


*  exp  [i(kz  +  k(r-z))]  ^ 

f  fmm«( 9  )  (1/kr)  exp[i (kr-r|&n  2kr)]  X 


m 

s 


(2.3-2) 

The  first  term  is  a  pure  Coulomb  wave  with  r|  =  ZpZte2/hk 

and  the  second  term  is  an  outgoing  spherical  wave  with  f  »  (0),  the 

scattering  amplitude  term  for  an  incoming  particle  scattered  through  a 

in  m  t 

centre  of  mass  angle  0;  while  X  and  X  are  the  spin-wave  functions 

s  s 

for  the  incoming  and  outgoing  particles. 


.A  partial  wave  expansion  for  Tm  in  the  Schrodinger  equation 

(H071)  with  total  anglar  momentum  J  and  orbital  angular  momentum  L  is 

written  asymptotically  as 

Y  ~  J  Att  (1/kr)  Rtt  Ftt  (2.3-3) 

m  r  +  00  L  JLm  JL  JLm 

where  the  spin-angle  functions  are  the  Fj^  terms  further  expressed  as 

r  .  m-m'  _» 

FJLm  =  ^  ^LS  ;m-m’  m'  I  Jm)YL(0,O)x®  (2.3-4) 

The  Clebsch-Gordan  (or  Wigner)  coefficient  terms, 

(LS;  m-ra'  m’|jm)  are  often  referred  to  as  vector-addition  coefficients. 


-  J 
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The  expansion  coefficient  terms }  are 

AjLm  =  iL  [47T (2L  +  1)]^  (LS;  Om|jm), 

(2.3-5) 

while  the  radial  wave  functions,  R  (r)  satisfy  the  equation 

J  L 

[d2/dr2  -  L(L  +  l)/r2  +  k2  -  2d/h2 
*  <  JL | U(r) | JL)>  ]  RJL(r)  =  0 

(2.3-6) 

where  U(r)  is  the  optical  potential  discussed  in  Section  2.1. 

The  solution  for  the  above  equation  is  given  in  terms  of  the 
Coulomb  functions  F  and  G  , 

Li  Li 

RjL (r )  ~  eiaL  [FL(t)  +ctL  +  iFL^r^] 

r  ■>  oo 

(2.3-7) 

and  the  absolute  Coulomb  phase  shifts,  defined  by  the  follow¬ 
ing  relation, 

exp(2iaL)  =  T(1  +  L  +  in)  /  T(1  +  L  -  in) 

(2.3-8) 

and  finally  the  nuclear  scattering  amplitudes,  further 
expressed  in  terms  of  the  nuclear  phase  shifts  6!^,  or 

=  (1/21)  exp  (216:5  -  1 


(2.3-9) 
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Asymptotically,  the  Coulomb  functions  and  are  defined 


by 


sin 


{  L  }  _ ►  {  }  (  kr-r|£n  2kr  -  Ltt/2  +aL  ) 

Gl  r  +  ° 


COS 


(2.3-10) 


Thus;  rewriting  equation  (2.3-3)  as  a  Coulomb  wave  function  and 
the  partial  wave  expansion  for  the  outgoing  sherical  wave, 


m  ^  -+  oo  m 


Coulomb 


+  Z  [4tt(2L  +  1)**  (LS;0m|jm) 
JL 


*(l/kr)  a  exp[i  (kr-Tl^n  2kr)] 

JL 

*  exp(2iCT^)  Z  (LS  ;  m-m  '  m  ’  [  Jm) 


..  T  1 

.  Y  m-m  ,  .  ..  m 

*  L  (  0,0)  x 


(2.3-11) 


tains 


and  rewriting  equation  (2.2.2)  in  a  similar  manner,  one  ob- 


m  t  00  1  m 


coulomb  +  £  (W<e)-fco)5mi?,)  (2>3_12) 

*  (l/kr)exp[i  (kr-T]£n  2  kr)]xm 


where  fc  is  the  Coulomb  scattering  amplitude 
fc(0)  =  -(o/2)csc2(0/2)  exp  (-2ir)£n  sin(  +  2iaQ) 


(2.3-13) 


Equating  co-efficients  from  equations  (2.3-10,  12),  one  obtains 
a  partial  wave  expansion  for  the  nuclear  scattering  amplitude 


f  ,  (0)  =  f  (0)  6  ,  +  E  [Att (2L+1) ] 

mm  c  mm  L  v  J 

J  Lj 


1/2 


*exp(2iGL)a^  (LS;  om|  Jm)  (LS  ;  m-m’m1  |  Jm)Y™  m  (0,0) 


(2.3-14) 


For  the  elastic  scattering  of  spin  ^  or  1  particles  from  a  spin  0 

nucleus,  the  expansions  for  f  .(0)  are  given  in  Appendix  A. 

mm 


The  differential  cross  section  for  incident  particles  is 
given  in  terms  of  the  scattering  amplitude  f  , (0) 

G(0)  =  [l/k2(2S+l)]  Tr  (f(0)f*(0)} 

2 


=  [l/k2 (2S+1) ]  Z  |  fm.m(6)| 


mm 


(2.3-15) 


Rewriting  a(0)  in  terms  of  the  scattering  amplitudes  f  , (0) 


mm 


found  in  Appendix  A  for  spins  ^  and  1  respectfully, 
S  =  h 

Q(0)  =  (l/k2)  (  |  f^l  2  +  |  2} 


(2 . 3-16a) 


S  =  1 


0(0)  =  (l/3k2)  {|fQ0| 2  +  2(|f1J2  +  |f1Q|  2  +  |fQ11 2  +  lf!_i|2) 


11 
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(2.3-16h) 
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CHAPTER  3 


EXPERIMENTAL 


3.1  Experimental  setup 

In  the  magpak  experimental  area  at  the  Nuclear  Research 
Centre,  the  incident  proton  or  deuteron  beam  from  the  Van  de  Graaff 
accelerator  was  collimated  by  two  adjustable  low-background  beam¬ 
defining  slits  upstream  and  an  antiscattering  slit  (Ro69)  located  at 
the  scattering  chamber,  (fig.  3.1)  Measurements  were  taken  in  an  18- 
inch  diameter  scattering  chamber,  where  inside  were  four  surface-barrier 
detectors,  that  is,  three  moveable  detectors  at  10  degree  intervals, 
located  at  a  distance  of  8  inches  from  the  centre  while  the  fourth 
detector,  9  inches  from  the  centre  and  stationary  at  30  degrees  to  the 
beam,  used  as  a  monitor.  An  important  feature  of  the  chamber  allowed 
the  detectors  to  be  moved  in  one  degree  steps  or  less.  In  front  of  each 
detector, circular  collimators  of  different  solid  angles  were  placed 
depending  on  the  count  rate.  By  observing  2.5  MeV  protons  elastically 
scattered  from  40Ca  the  relative  solid  angles  of  the  collimators  were 
measured.  The  output  analog  pulses  from  each  of  the  detectors  went 
through  a  preamp-amp-ADC-DDP516  configuration  for  real-time  data 
acquisition  and  finally  into  a  XDS940  computer  for  off-line  data  analysis. 

The  targets  were  made  by  evaporating  approximately  20  yg/cm2 
of  natural  calcium  metal  (96.6%  40Ca)  onto  a  carbon-formvar  backing.  The 
thin  layer  of  carbon  was  used  to  dissipate  heat  from  the  beam  while 
the  formvar  was  used  for  mechanical  strength  of  the  target.  Because  of 
the  chemical  properties  of  calcium,  that  is,  oxidation  of  the  metal, 
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Schematic  of  Experimental  Arrangement 


the  problem  of  eliminating  oxygen  completely  was  impossible 
with  the  present  evaporating  sys tern-target  chamber  method  since 
air  was  always  present  at  some  stage  of  the  preparation. 
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3.2  Relative  cross-sections 

Cross-sectional  data  was  measured  in  five  degree  steps  or 
less,  in  the  angular  distribution  range,  7-160  degrees,  and  many  of  the 
experimental  runs  were  repeated  to  obtain  statistical  results  with  the 
following  corrections: 

i)  background  and  multiple  scattering 

The  probability  that  n  events  will  occur  in  any  one 
channel  of  an  ADC  energy  spectrum  is  a  Poisson  distribution.  If  the 
number  of  events  or  counts  under  a  typical  elastic  peak  in  the  energy 
spectrum  is  large,  then  the  approximation  to  a  gaussian  function, 

P(x)  =  (A/  2a)exp[-(%)  ((x-tO/g)2] 

(3.2-1) 

with  parameters  A  for  the  area,  O  for  the  standard  deviation  in  A  and 
x  is  the  c.m.  energy  centroid;  is  a  standard  least  squares  peak  fitting 
procedure.  The  independent  weights  w^  assigned  in  the  equation  below 
for  the  minimizing  function  A, 

A  =  Z  w.  [PTH(x)-PGXp(x)] 

i  1  1  (3.2-2) 

are  the  square  root  of  the  number  of  counts  for  each  channel.  The  peak 
fitting  program  has  the  following  features: 

a)  background  subtraction 

b)  separation  of  12C  and  160  peaks  from  40Ca 

c)  separation  of  multiple  scattering  events  (i.e.  in  the 
forward  angle  region)  which  produce  a  skew  gaussian 

d)  evaluation  of  a  physical  statistical  error  for  the 


area  under  the  peak. 


♦ 
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'  ii)  dead-time  corrections  for  the  ADC's. 

If  the  logic  gate  is  open,  an  ADC  will  accept  a  random 
pulse  from  a  detector  and  become  busy,  i.e.  will  not  accept  any  more 
pulses.  The  analog  pulse  is  converted  to  a  digital  pulse.  The  gate 


is  closed  and  the  ADC 

is  ready  for  the  computer.  When  the  information 

is  dumped,  the  ADC  is 

reset  and  the  gate  is  open  to  accept  another 

random  pulse. 

Since  the  count  rate  varies  approximately  as  l/sin4 (0 /2)  , 
the  dead  time  for  the  forward  angles  may  be  greater  than  at  the 
backward  angles  where-  the  count  rate  is  not  large. 

The  relative  peak  sum  was  normalized  to  the  monitor  located 
30  degrees  to  the  beam  with  the  following  correction. 

Psum  =  Mav  Dsum/DLT 

Msum/MLT  (3.2-3) 

where  the  following  symbols  are  defined  below: 


Psum  - 

relative  peak  sum, 

Dsum  - 

detector  peak  sum, 

Ms  urn  - 

monitor  peak  sum. 

DLT  - 

detector  livetime, 

MLT  - 

monitor  livetime, 

Mav 

average  monitor  sum. 

* 
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iii)  Solid  angle  corrections 

. 

Since  the  count  rate  varies  approximately  as  l/sin4(0/2), 
collimators  of  different  solid  angle  were  placed  in  front  of  each 
detector.  Normalization  was  done  using  2.5  MeV  protons  elastically 
scattered  from  40Ca  at  30  degrees. 


iv)  target  thickness,  including  deterioration  and 
beam  stability  using  low  energy  protons. 

v)  target-angle  corrections  at  90  degrees. 


The  target  was  at  330  degrees  to  the  beam  for  the  forward 
angles  and  changed  to  30  degrees  for  the  backward  angles. 


vi)  conversion  of  the  lab  system  to  the  centre  of 
mass  system 


da(Q) 

dft 


cm 


m  mUm/m  )  2  -  s±nZy¥  da  O')  i 

p  t  t  p  _  I 

(m  +  m  )2 (E, /E.  )  d ^  Lab 

^  p  t  i  me 

(3.2-4) 


where  ip  is  the  laboratory  angle  and  E.  is  the 

inc 

incident  particle  energy,  while  is  the  energy 
of  the  particle  scattered  from  a  target 
nucleus.  The  other  parameters  are  defined  in 
Section  2.1. 

vii)  since  the  experimental  runs  were  short,  i.e.  5-8 
minutes,  many  of  the  runs  were  repeated  to  obtain 
good  statistics.  Defining  x  as  the  mean  peak  sum 


■ 


- 


area, 


x 


E  w.x./E  w.  ,  i=l,N 
11  1 
i  i 


(3.2-4) 


and  the  weight  w. 


is  written  as 


w.  =  1/  A  x? 

l  l 


(3.2-5) 


is  related  to  the  error  in  the  peak  sum. 
the  definition  of  chi-squared 


£  w±(x  -  x±) 2  , 


(3.2-6) 


therefore  the  error  in  Ax  is 

&x  =  {(1/Ewp  (x2/(n-l))}1/2 

(3.2-7) 


From 


where  n  is  the  number  of  degrees  of  freedom. 
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3.3  Absolute  Cross  Sections 

From  the  elastic  scattering  of  low-energy  protons  in  the 
2.5  -  3.5  MeV  range,  the  absolute  cross  sections  should  be  nearly 
Rutherford  in  the  forward -angle  region  and  the  relative-to-absolute 
conversion  factor  should  not  be  difficult  to  obtain  for  a  target  nucleus  such 
as  calcium.  However,  for  the  light  nuclei,  carbon  and  oxygen,  nuclear 
contributions  (BL65,  HA62  and  RE56)  provide  difficulties  in  finding 
satisfactory  conversion  factors;  as  a  result, the  absolute  cross-section 
error  bars  are  large.  To  show  the  above  behavior,  theoretical  curves 
are  generated  (figs.  3.1  -  3.2)  using  an  optical-model  search  program  on 
the  proton  experimental  data  starting  with  the  initial  parameters  in 
Table  3.1  below. 


Table  3.1  Initial  Parameters  for  12C  and  160  Proton  Data 


V 

V 

Ao 

WD 

Rw 

vso 

Rso 

Aso 

(MeV) 

(fm) 

(fm) 

(MeV) 

(fm) 

(fm) 

(MeV) 

(fm) 

(fm) 

50.0 

1.25 

0.65 

8.5 

1.25 

0.52 

9.0 

1.25 

0.64 
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Searching  only  on  the  potentials  W  and  V  n  with  the  computer 

U  DU 

* 

program  SNOOPY  gave  a  nonphysical  range  of  values  for  both.  W^  was  be¬ 
tween  7.21  -  17.23  MeV  and  the  range  of  the  spin  orbit  potential 
was  from  1.75  -  35.28  MeV.  The  geometrical  potential  parameters  were 
calculated  and  gave  no  improvement  to  the  fit.  The  nonphysical 
potential  parameters  calculated  by  the  search  program  SNOOPY  were 
due  to 


i)  few  experimental  data  points, all  of  which  were  in  the 
forward  angle-region; 

*  i  a  i  r 

ii)  the  resonance  structure  present  for  both  C  and  0 
nuclei  (BL65,  HA62  and  RE56)  in  the  low-energy  proton 
spectrum  region;  and  finally 

iii)  the  optical  model  breaks  down  below  5.0  MeV. 

The  only  reasonable  parameters  to  use  for  calculations  of  the 
optical  model  are  found  in  Table  3.1,  obtained  from  references 
H071  and  PE63.  From  figures  3.1  and  3.2,  one  finds  the  relative- 
to-absolute  conversion  factors  are  difficult  to  acquire  with  any 
degree  of  accuracy;  therefore,  the  absolute  cross  sections  for  C 
and  160  experimental  data  were  obtained  from  normalizations  to  5.0  MeV 
deuteron  data  found  elsewhere  in  the  literature  (DA70  and  OH63) . 

/ 

*  Obtained  from  W.  Haeberli,  University  of  Wisconsin,  Madison, 
Wisconsin . 


Differential  cross-section  as  ratio  to  Rutherford  cross-section 
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Centre -of -mass  angle 

Figure  3.1  2.5,  3.0  and  3.5  MeV  12C(p,p)12C 


Differential  cross-section  as  to  ratio  to  Rutherford  cross-section 
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Cen tre-of-raass  angle 

Figure  3.2.  2.3,  3.0,  and  3.5  MeV  160(p,p)160 


CHAPTER  4 


RESULTS  AND  DISCUSSION 

A.l  Op tical- model  Analysis  of  5.0  MeV  Data. 

From  the  elastic  cross-section  measurements,  an  attempt  to  find 
the  optical  potentials  that  fit  the  experimental  data  was  made  using  a 
computer  program  called  SNOOPY.  The  fitting  procedure  is  as  follows: 

i)  optimize  fits  to  the  experimental  data  starting  first  with  the 

potentials  V,  W^,  and  while  leaving  the  geometrical  parameters  fixedj 

ii)  use  the  maximum  number  of  partial  waves  so  the  nuclear  contributions 
from  all  L  values  are  taken  into  account ; 

iii)  vary  "  the  geometrical' parameters  R,.,  A  ,  R  ,  A  ,  R  and  A 

U  U  W  W  bU  bU 

each  separately  leaving  the  rest  fixed  ; 

iv)  finally  simultaneous  vary  all  the  parameters  together. 

From  i)  and  iii),  the  general  fitting  behavior  of  the  parameters 
may  be  studied,  i.e.  fits  the  backward  angles  and  V  SO  fits  the  dip 
around  90  degrees.  Searching  on  all  the  paramets  in  iv)  at  once  should 
give  little  change  in  both  the  parameter  and  chi-squared  values.  Using 
only  a  few  L-values  in  ii) ,  say  computed  from  the  classical  result 
L  =  r  x  one  finds  sometimes  the  parameters  varied  gave  nonphysical 
results . 

The  parameters  may  not  be  physically  correct  even  though  they 
fit  the  cross-sectional  data.  For  example,  if  there  is  a  resonance 
structure  present  one  would  expect  rapid  variation  of  the  optical-model 
parameters  with  energy.  Therefore  the  proper  choice  is  to  use  parameters 
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changing  smoothly  with  energy.  One  should  also  remember  that  the  optical 
model  description  for  elastic  scattering  gives  only  the  asymptotic 
behavior  of  the  wave  functions,  but  indicates  little  about  what  goes  on 
inside  the  nucleus. 

The  initial  parameters  searched  for  the  elastic  scattering  of 
5.0  MeV  deuterons  on  12C,  160  and  l4°Ca  nuclei  are  shown  in  Table  4.1 
below. 


Table  4.1 


INITIAL  OPTICAL -MODEL  PARAMETERS 


V 

(MeV) 


100 


Ro 

(fm) 


1.05 


0 

(fm) 


0.85 


WD  \ 

(MeV) (fm) 


8.5 


1. 66 


\  vso  R~  A 


(fm) 


0.53 


(MeV ) 


9.0 


SO 

(fm) 


0.9 


SO 

(fm) 


0.65 


The  function  y2 ,  defined  as 

i  N 

x2  =  (  /N)£  [(da(0)dfi)TH  -  (da(e)/dfi)EXp)/(Ada(e)/dQ)]2 

i=l 

was  used  to  compare  the  fits  in  terms  of  the  optical-model  parameters. 

Here,  (da(0)/dS7)  and  (da(0)/dfi)  are  the  calculated  and  experimen- 
tal  absolute  cross  sections,  respectively,  for  the  centre-of~mass 
angle  0.  while  the  weight  assigned,  (Ada(0)/df2)  is  the  experimen- 

tal  error.  The  above  expression  for  y2  assumes  the  errors  are  uncorre¬ 
lated  for  each  data  point.  However,  errors  from  the  elastic  scattering  events 


' 


2  4 


come  not  only  from  the  statistics  of  the  number  of  points  or  area 
in  the  elastic  peak,  or  the  number  of  repeated  runs,  but  also 
from  measurements  of  target  thickness,  solid  angle,  dead- time 
corrections  and  so  forth  which  may  or  may  not  be  correlated  with 
angle  0. 

The  present  analysis  uses  an  automatic  parameter  search  routine 
which  minimizes  chi-square  with  respect  to  any  desired  number  of 
parameters.  The  numerical  value  of  X2  obtained  does  not  give  a  very 
definite  idea  of  the  fit  while  the  use  of  X  as  a  measure  of  goodness 
of  fit  is  more  meaningful.  Frequently, the  automatic  search  ends  in  a 

.  ,  o 

secondary  minimum  or  valley  of  the  X  surface  in  parameter  space, 
particularly  when  the  number  of  parameters  varied  is  large  or  when  only 
a  few  partial  waves  are  used. 


. 


' 
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4.2  5.0  MeV  12C(d,d)12C  and  160(d>d)160. 


The  theoretical  curves  in  figures  4.1  -  4.2  are  generated  from 
the  optical  model  parameters  shown  in  Tables  4-2  and  4.3  below. 


Table  4.2  OPTICAL  MODEL  PARAMETERS  FOR  5.0  MeV  12C(d,d)12C 


V 

Ro 

Ao 

WD 

*w 

*w 

vso 

Rso 

Aso 

(MeV) 

(fm) 

(fm) 

(MeV) 

(fm) 

(fm) 

(MeV) 

(fm) 

(fm) 

104.96 

1.136 

0.85 

5.39 

1.612 

0.53 

8.89 

0.792 

0.64 

Table  4.3  OPTICAL  MODEL  PARAMETERS  FOR  5.0  MeV  160(d,d)160 


V 

Ro 

Ao 

WD 

% 

AW 

vso 

Rso 

Aso 

(MeV) 

(fm) 

(fm) 

(MeV) 

(fm) 

(fm) 

(MeV) 

(fm) 

(fm) 

102.35 

1.046 

0.85 

4.81 

1.992 

0.53 

11.99 

0.70 

0.64 

The  theoretical  optical  dip  about  20  degrees  found  in  figures 
4.1  -  4.2  is  important  since  the  deviation  from  Rutherford  scattering 
for  both  12C  and  160  nuclei  at  5.0  MeV  is  shown  here.  For  all  sets  of 
parameters  searched  on,  even  varying  the  number  of  partial  waves  one 
finds  that  the  nuclear  part  of  the  scattering  amplitude  dominates;  in  other 
words,  the  dip  was  always  present.  Since, from  the  energy  spectrum,  the 
centre  of  mass  of  both  12C  and  1G0,  elastic  peaks  shifts  into  40Ca,  it  becomes 
rather  difficult  to  analyse  for  the  nuclear  contributions  experimentally. 


•: 


Differential  cross-section  as  ratio  to  Rutherford  cross-section 


20  40  60  80  100  120  14C  160  180 

Cen tre-of-raass  angle 
Figure  4.1  5.0  MeV  12C(d,d)12C 


Differential  cross-section  as  ratio  to  Rutherford  cross-section 
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Cen tre-of-mass  angle 

Figure  4.2  5.0  MeV  160(d,d)160 
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Even  if  one  could  predict  the  nuclear  part  of  the  scattering  amplitude, 
there  would  still  be  the  interference  term  between  the  nuclear  and 
Coulomb  parts  shown  below 


+ 


— i  (4> 

+  e 


N 


H  |fj 


2 


=  lfNl2  +  2lfJ  lfJ  cos^ 


N 


N 


-  *c) 


where 

N 


is  the  interference  angle. 
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4.3  5.0  MeV  40Ca  (d,d) 40Ca. 

The  results  from  the  optical  model  calculations  are  shown  in 
fig.  4.3  and  in  the  table  below. 

Table  4.4  OPTICAL  MODEL  PARAMETERS  FOR  5.0  MeV  4 °Ca (d, d) 4 °Ca 


V 

(MeV) 

Ro 

(fm) 

Ao 

(fm) 

md 

(MeV) 

(fm) 

(fm) 

vso 

(MeV) 

Rso 

(fm) 

Aso 

(fm) 

103.35 

1.053 

0.934 

9.82 

0.655 

0.511 

8.5 

1.302 

0.725 

The  sensitivity  of  the  optical  model  parameters  particularly 
to  small  deviations  from  Rutherford  around  15  degrees  in  figure  4.3  is 
important.  The  physical  potential  parameters  obtained  from  the  optical 
model  calculations  were  not  sensitive  in  the  above  forward  angle 
region.  Using  only  the  experimental  data  points  less  than  90  degrees 
again  predicted  the  same  behavior. 


’  pH|  '  '  .  v-'  1  ‘  1 1  '8  } 


Differential  cross-section  as  ratio  to  Rutherford  cross-section 
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Centre-of-mass  angle 


Figure  4.3  5.0  MeV 


40Ca(d,d) 4  °Ca 
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4.4  Summary  and  Conclusions 

From  the  analysis  of  the  present  data  it  may  be  concluded  that  the 
elastic  scattering  of  5.0  MeV  deuterons  on  40Ca  is  not  pure  Rutherford 
scattering  in  behavior  (nuclear  bump  around  15  degrees) .  Also  solving  the 
problem  of  separating  the  other  elastic  peaks  from  one  another  in  the 
forward  angle  region  represents  a  problem  unsolved. 

Two  proposals  may  be  suggested: 
i)  Run  5.0  MeV  deuterons  on  a  carbon  backing  at  the  forward - 

angular  region.  From  the  elastic  scattering  of  2.0  MeV  protons 
on  C,  the  approximate  amount  of  0  may  be  found.  For  some 
angles  where  the  elastic  peaks  are  separate  compare  with  the 
angular  distribution  results  of  40Ca  on  a  carbon  backing  for 
oxygen  contamination. 

ii)  Make  an  oxygen  gas  cell  for  the  target  chamber  and  obtain 
an  angular  distribution  from  2.0  MeV  protons  and  5.0  MeV 


deuterons . 


.. 
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APPENDIX  I 


Nuclear  Scattering  Amplitudes  for 
Spin  ^  and  1  Particles 


From  Section  2.3,  a  partial  wave  expansion  of 

the  scattering  amplitude  was  obtained  for  the  f  t (0)  term 

mm 


in  equation  2.3-14.  The  explicit  expressions  for  f 


(0) 


are  now  given  below. 


i)  S 


-  L 


f11(0)=fC0)+Ee  L  [  (L+l)  +  La^  ^]PT  (cos  0) 

c  L  L  L 

Li 


2  ia 


fv_k(0)  = 

-"2 - 2 


Ee 

L 


L  r  Li  +  ^2  L  —  %  1 


[a 


-  aT  ]  PT  (cos  0) 
L  Li 


ii)  S  =  1 


foo(e) 


2 ia  L  L+l  L— 1  . 

f  (0)  +  Ee  [(I/t-1)  a  +  La  ]P.  (cos  0) 

Q  Li  Lj  -Li 

L 


fllce) 


f  (0)  +  Ee 

C  L 


2iC 


L  ih)  [  (L+2)a^+1  +  ( 2L+ 1 ) a 

Li 


L 

L 


+  (L-l)a^  1]P  (cos  0) 

Li  Li 


2  ia  .  , 

f  (6)  =  Ee  L(l//2)  [ct^  -  a  ]  P  (cos  0) 

10  L  Li  Li 

Li 
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34 


2  ia 


foiC6>  -  Ee 

L 


[1//2  L(L+1)  J  I-LCL+2)  a 


L+l 


+  (2L+l)a^  +  (L+l)  (L-l) 1]  p 1  (cos  0) 
L  L  l 

2  ia 

f,  -j  C  0 )  =  Ee  L  [  %L(L+1)  ]  [La^  -  (2L+l)a^ 

-L  ~  J-  L 

L 


L- 1  2 

+  (L+l )  a  ]  P  ( cos  0 ) 
L  L 


